A steady stratified purely azimuthal flow representing the Antarctic Circumpolar Current

We construct an explicit steady stratified purely azimuthal flow for the governing equations of geophysical fluid dynamics. These equations are considered in a setting that applies to the Antarctic Circumpolar Current, accounting for eddy viscosity and forcing terms.


Introduction
Thorough analytical investigations of exact solutions to the fully nonlinear governing equations of geophysical fluid dynamics (GFD) represent an extensive and active research area, which was initiated by Constantin [1][2][3][4] and Constantin and Johnson [5][6][7][8]. Following this approach, we construct a unidirectional flow satisfying GFD considered in the so-called f -plane at the 45th parallel south, enhanced with an eddy viscosity term and a forcing term, and equipped with appropriate boundary conditions. We propose this specific flow for representing the gross dynamics of the Antarctic Circumpolar Current (ACC)-the World's longest and strongest Ocean current.
ACC has no continental barriers: it encircles Antarctica along a 23 000 km path around the polar axis towards East at latitudes between 40°to 60°, see Fig. 1. It thereby links the Atlantic, Pacific and Indian Oceans making it the most important oceanic current in the Earth's climate system. ACC has a rich and complicated structure. Many factors contribute to its complex behavior-the most important driver are the strong westerly winds in the Southern Ocean region. In addition to that there exist mesoscale eddies of a size up to 100 km, which transport the wind-induced surface stress to the bottom and also enable meridional mass transport; there are sharp changes in water density due to variations in temperature and salinity-known as fronts or jets-located at ACC's boundaries (see Fig. 1); ACC is strongly constrained by the bottom topography; there are observed variations in time such as the Antarctic Circumpolar Wave; etc. We refer to [9,11,12,16,18,19,24,25,27,28] for further information about the geophysical aspects and modeling as well as observational data and simulations for ACC.
From an analytical perspective one is forced to largely, yet reasonably, simplify the geophysical scenario to obtain a tractable model, which-in the ideal case-exhibits exact and explicit solutions opening the path for an in-depth analysis. Thus we do not account for all of the before mentioned phenomena, but assume a steady flow in purely azimuthal direction, which is vertically bounded by a flat bottom and a flat ocean surface. By considering Euler's equation of motion in the f -plane, we obtain a valid approximation of the Coriolis effects close to the 45th parallel south; in this way the Earth's curvature is neglected and no boundaries in the meridional direction are assumed. To account for the transportation effects of mesoscale eddies we equip the system with an eddy viscosity term; furthermore we include a forcing term to ensure the dynamical balance of the flow. Both pressure and wind stress are prescribed on the ocean surface; a no-slip boundary condition is assumed for the ocean bed.
A similar setting has recently been considered in [22], where an explicit solution in terms of a given viscosity function was presented. In this note at hand we extend these results to stratified flows, i.e. we do account for variations of the water density (with depth and latitude). The established explicit solution is an analytic function of both the viscosity function and the density distribution. A collection of numerous recent related analytical studies can be found in [7,13,14,17,20,22] and the references therein.

Model under study
We introduce the governing equations for geophysical ocean flows that set the basis for our study: we take into account the effects of the Earth's rotation on the fluid body by choosing a rotating framework with the origin at a point on the Earth's surface. Accordingly, we will work with Cartesian coordinates (x, y, z), where x denotes the direction of increasing longitude, y is the direction of increasing latitude and z represents the local vertical, respectively. Denoting with t the time variable, and with (u(x, y, z, t), v(x, y, z, t), w(x, y, z, t)) the velocity field, the governing equations for inviscid and incompressible geophysical ocean flows at latitude φ are  where u = (u, v, w). Here, P = P (x, y, z, t) denotes the pressure field, ρ = ρ(y, z) serves a (prescribed) density distribution to account for density gradients, mainly caused by variations of salinity and temperature, Ω ≈ 7.29 · 10 −5 rad s −1 is the Earth's (constant) rotational speed around the polar axis toward the East, and g ≈ 9.82 m/s 2 denotes the gravitational constant. The governing equations (2.1)-(2.2) hold throughout the fluid domain, which lies between the rigid flat bed at z = −d (d is the constant water depth) and above the flat surface situated at z = 0. Our aim is to derive exact formulas for purely azimuthal flows (i.e. v = w = 0) in the region of the ACC (i.e. we take φ = −π/4, which corresponds to the 45th parallel south). Furthermore we incorporate the transfer of the wind-generated surface stress to the bottom, which is due to the presence of mesoscale eddies, by adding a viscosity term of the form (νu z ) z to the right hand side of the first equation in (2.1); the coefficient ν = ν(z) is a smooth function of depth being strictly greater than some positive constant, see [16]. The classical model of uniform eddy viscosity is due to [23]. We follow the more realistic approach with a depth dependent viscosity function as it was introduced in [10]. Finally, we include a forcing term F = F (y, z) to guarantee non-trivial solutions, cf. [7,15,22].
According to the previous considerations, we will consider the following set of equations governing the ACC as a geophysical flow model in the purely azimuthal direction: where τ (y, z) := ρνu z represents the wind stress; i.e. we assume a no-slip bottom and constant pressure as well as wind stress at the surface. Proof. Utilizing (2.6) we differentiate by x in (2.3)-(2.5) and obtain that

Explicit solution
From the condition P = P atm on the surface z = 0 we obtain that P x = 0 on z = 0. Thus, from (3.3) we see that P x = 0 in D. Hence, equation (2.3) becomes (νu z ) z = 0, which implies (making also use of (2.6)) that there is a funtion y → B(y) such that The latter equation and condition (2.9) imply that ν(0)u z (y, 0) = τ 0 ρ(y, 0) = B(y) for all y.
Integrating with respect to z in (2.5) we obtain that P (y, z) = √ 2Ω Remark 3.2. One immediate consequence of (3.1) is that the vorticity vector associated with the flow (3.1) is (0, u z , −u y ) has a non-vanishing second and third component. This represents a marked difference, if compared with the case of homogeneous flows (considered in [22]) where only the middle component u z survives, the first and the third being zero because of the lack of y dependence of ρ. Thus, allowing for significant variations in density leads to solutions that exhibit substantial shear not only in the vertical direction but also in the latitudinal direction, as well.